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ABSTRACT 

We  prove  a  cell  entropy  inequality  for  a  class  of  high  order  discontinuous  Galerkin  finite 
element  methods  approximating  conservation  laws,  which  implies  convergence  for  the  one 
dimensional  scalar  convex  case. 
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1  Introduction 


Iu  [3]  Cockburu  and  Shu  defined  a  class  of  discontinuous  Galerkin  finite  element  methods 
for  conservation  laws  (the  multi-dimensional  case  was  discussed  in  [4]): 

tti  +  div /(u)  =  0  (1.1) 

The  entropy  solution  of  (1.1)  also  satisfies 

U(u)t  +  div  F(u)  <  0  (1.2) 

in  the  distribution  sense,  for  any  convex  function  U (u)  and  consistent  entropy  flux  F(u)  sat¬ 
isfying  F'(u)  =  U'(u)f'(u).  The  scheme  is  obtained  in  the  following  way:  first  a  triangulation 
Th  is  chosen.  In  one  space  dimension  Th  is  just  a  collection  of  subintervals  Ij  =  (Xj_i,  xJ+i), 
which  are  not  necessarily  of  the  same  length.  We  will  use  the  notation  Aij  =  xJ+i  —  £j_i, 
h  =  maxj  Aij  and  xj  =  |(xJ_i  -f  xJ+i).  In  two  and  more  space  dimensions  7*  is  a  collection 
of  triangles,  tetrahedrons  or  other  simple  geometric  objects.  The  solution  space  Vh  is  defined 
as  the  collection  of  all  piecewise  polynomials  of  degree  up  to  r  for  a  (r  -f  1  )-th  order  method. 
The  functions  in  Vh  are  allowed  to  have  discontinuities  across  element  interface.  The  conser¬ 
vation  law  (1.1)  is  multiplied  by  a  test  function  v  €  Vh,  integrated  over  an  element  K  €  Th, 
and  formally  integrated  by  parts  to  shift  the  spatial  derivatives  from  /(u)  to  v.  The  result 
is: 


f  ut(x,t)v(x)dx  +  ^2  j  f(u(x,t))-nv(x)ds  —  f  f(u(x,t))  •  grad  v(x)dx  =  0.  (1.3) 

JK  eedKJe  Jk 

where  dK  is  the  boundary  of  the  element  K  and  n  is  the  unit  outward  normal  vector.  This 
procedure  up  to  now  is  just  the  standard  discontinuous  Galerkin  method.  Two  features  from 
the  finite  difference  methodology  are  then  used  in  [3],  [4]:  one  is  the  monotone  or  E  flux 
(approximate  Riemann  solver  for  systems)  to  define  the  trace  of  the  physical  flux  /(u)  at 
the  cell  interface;  the  other  is  a  local  nonlinear  limiter  which  limits  the  values  of  u  at  the 
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cell  interface  by  the  differences  of  the  means  of  u  over  elements.  With  these  two  ingredients, 
it  was  proven  in  [3]  and  [4]  that  the  schemes  are  formally  high  order  accurate,  measured  by 
local  truncation  errors  in  smooth  regions  including  at  extrema,  are  total  variation  bounded 
for  one  space  dimension,  and  are  maximum  norm  bounded  for  any  space  dimensions.  The 
only  thing  missing  is  the  entropy  condition:  in  [3],  we  were  able  to  prove  entropy  consistency 
for  the  square  entropy  U(u)  =  for  one  dimensional  convex  f(u)  with  a  ^-independent 
modification  to  the  scheme,  following  the  idea  of  Osher  [12].  We  were  also  able  to  prove 
entropy  consistency  for  all  convex  entropies  with  a  /i-dependent  modification  to  the  scheme. 
The  /j-dependent  limiters  make  the  proof  of  high  order  schemes  easy,  but  it  is  not  very  desir¬ 
able  for  practical  computations,  because  it  usually  limits  the  slope  near  the  discontinuities 
stronger  than  necessary,  and  essentially  flattens  the  solution  there  to  piecewise  constants  for 
fine  mesh.  It  also  destroys  the  self-similarity  of  the  scheme. 

Entropy  condition  seems  difficult  to  prove  for  high  order  finite  difference  schemes.  Osher 
and  Tadmor  proved  [13]  that  finite  difference  schemes  (which  evolve  only  the  mean)  which 
satisfy  cell  entropy  inequalities  for  all  convex  entropies  can  be  at  most  first  order  accurate. 
Even  for  one  entropy  inequality  (say  for  the  square  entropy),  for  one  space  dimension  and 
for  convex  /(u),  the  proof  is  extremely  elusive  if  one  does  not  modify  the  scheme.  Osher  [11] 
aud  Nessyahu  and  Tadmor  [10]  were  able  to  get  such  cell  entropy  inequalities,  for  the  square 
entropy  and  for  the  second  order  MUSCL  scheme,  with  some  ^-independent  modifications; 
Yang  [14]  was  able  to  prove  convergence  of  unmodulated  second  order  MUSCL  scheme 
using  a  global  analysis  rather  than  relying  solely  on  cell  entropy  inequalities;  Lions  aud 
Souganidis  [8]  proved  convergence  of  second  order  MUSCL  scheme  for  steady  state  Hamilton- 
Jacobi  equations  and  conservation  laws.  There  are  also  many  results  which  prove  entropy 
consistency  and/or  convergence  using  Ji-dependent  limiters  or  modifications  for  high  order 
schemes.  For  example  Coquel  and  LeFloch  [7]  for  finite  difference;  Johnson,  Szepessy  and 
Hansbo  [9]  for  streamline  diffusion  finite  element;  Cockburn,  Coquel  and  LeFloch  [5]  and 
Cockburn  and  Gremaud  [6]  for  high  order  finite  volume,  etrcemline  diffusion  or  discontinuous 
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Galerkin  schemes  with  /i-dependent  “shock  capturing”  terms.  These  results  are  usually  more 
general  (multi  space  dimensions,  noil-convex  fluxes,  etc.).  However,  as  we  have  indicated 
before,  ^-dependent  limiters  or  modifications  should  be  avoided  for  practical  computations  if 
possible.  More  recently,  Bouchut  and  Perthame  [2]  obtained  a  second  order  one  dimensional 
scheme  which  is  consistent  with  all  entropy  conditions  and  does  not  use  ^-dependent  limiter. 
A  key  ingredient  of  [2]  is  to  evolve  both  the  mean  and  the  slope  and  use  the  whole  function 
(not  just  the  mean)  to  obtain  cell  entropy  inequalities.  Discontinuous  Galerkin  methods  also 
fall  into  this  category  (evolving  the  whole  polynomial  in  a  cell  rather  than  just  the  mean), 
and  we  use  this  fact  strongly  in  this  paper. 

In  Section  2  we  prove  a  cell  entropy  inequality  using  the  square  entropy  U(u)  =  y 
for  the  unmodulated  semi-discrete  discontinuous  Galerkin  method  of  [3],  [4],  The  proof  is 
remarkably  simple  and  does  not  even  use  any  nonlinear  limiters.  It  works  thus  for  any  spatial 
order  of  accuracy.  As  far  as  we  know,  all  the  previous  cell  entropy  inequalities  without  h- 
dependent  limiters  must  restrict  the  slope  of  a  function  by  -g^-minmod(A+Uj,  A _Uj)  (which 
means  that  the  difference  at  the  interface,  u++1  —  u~+  L,  must  be  of  the  same  sign  as  that  of 
iij+i  —  uj,  i.e.,  no  “sawtooth”  is  allowed  in  the  reconstruction),  hence  cannot  be  higher  than 
second  order  accurate.  The  result  in  this  paper  illustrates  the  potential  of  discontinuous 
Galerkin  methods,  or  equivalent  Hermite  type  finite  difference/finite  volume  type  methods, 
which  evolve  the  whole  polynomial  in  the  cell  rather  than  just  the  mean.  Time  discretization 
is  discussed  in  Section  3. 

2  Cell  Entropy  Inequality  for  the  Square  Entropy 

The  discontinuous  Galerkin  scheme  in  one  space  dimension,  defined  in  [3],  in  its  semi-discrete 
form  without  slope  limiting,  is  the  following:  Find  u(-,t)  E  14,  such  that,  for  all  v  €  14  and 
all  subinterval  ly 


JI}  ««(*.  <)»(*)<**  +  hJ+i(t)v(x-+J  -  hj_  i(t)v(x+_L)  -  Ji  f(u(x,t))vx(x)dx  =  0  (2.1) 
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Here,  hj+^(t)  =  h  i(xj+i,<),u(xt+i,f)^  is  a  Lipschitz  continuous  monotone  flux  (i.e.  h  is 
non-decreasing  in  the  first  argument  and  non-increasing  in  the  second  argument),  or  more 
generally  an  .E-flux  as  defined  by  Osher  [11]: 

(/i(u-,u+)  —  /(«))  (u+  —  u_)  <  0  (2.2) 

for  all  u  between  u~  and  u+.  Some  examples  of  the  commonly  used  monotone  fluxes  can  be 
found  in,  e.g.,  [3]. 

If  we  take  v(x)  =  u(x,t)  in  (2.1),  we  get 


i,,(- '¥) 

We  define 


dx  + h  ,i(f)it(x ..t,t)-h ,_i(t)u(xt  ,,<)-  f  f(u(x,t))ux(x,t)dx  =  0  (2.3) 

J  a  Jtj  j  a  j  j  Jij 


g(u)  =  J  f{u)du 

and  rewrite  (2.3)  as 


(2.4) 


Jh  *  +  Fi+i(t)  -  FH(t)  +  Aj(l)  =  0  (2.5) 

where 

Ei+i(f)  =  hj+i(t)u(x'+±,t)  -  g(u(x-+±,t))  (2.6) 

is  consistent  with  the  entropy  flux  for  the  square  entropy: 


F(u)  =  J  f'(u)udu  =  f(u)u  -  J  f(u)du  =  f(u)u  -  g(u), 


(2.7) 


and 


Aj{t)  =  -h (u(x+_±,t)  -  u(x'_L,t)j  +y(u(*t_i,<))  -<7(u(xj_i,t)) 
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(V*(0  - /(O)  ^  0 


(2.8) 


where  we  have  used  the  mean  value  theorem  and  the  definition  (2.4)  of  <7(u)  in  the  second 

equality,  £  is  between  u(x~  k ,  t)  and  u(x+  i ,  t ),  and  the  last  inequality  is  due  to  the  property 

1  1  3  2 

of  E-flux  (2.2). 

We  have  thus  proven  the  cell  entropy  inequality 


u2(x,  t) 
.  2 


dx  +  Fj+Ut)~  <  0 


for  the  square  entropy  U(u)  =  Notice  that  we  do  not  need  any  nonlinear  limiting  at 
this  stage.  However,  nonlinear  limiting  as  introduced  in  [3]  and  [4]  will  not  destroy  this  cell 
entropy  inequality  (see  next  section).  The  cell  entropy  inequality  (2.9)  trivially  implies  L2 


stability  of  the  scheme  (again  without  even  using  the  nonlinear  limiting): 


(2.10) 


but  it  is  much  stronger.  For  example,  if  f(u)  is  convex  and  we  use  the  nonlinear  limiting  [3] 
to  obtain  a  total  variation  boundedness  for  the  solution,  we  will  have  convergence  towards 
the  unique  entropy  solution. 

The  same  entropy  inequality  can  be  obtained  for  multi-space  dimensions  with  arbitrary 
triangulations: 


/  <>*+  E  f  F,.K(x,t)dS<  0  (2.11) 

JK  \  1  It  eedKJ' 

where  Fe,K  is  consistent  with  F  •  n#  for  the  entropy  flux  F  in  (1.2)  and  the  outward  normal 
uk,  and  Fe,K  =  —Fe<K>  for  the  two  neighboring  elements  Kf]K'  =  {e}  (conservation).  We 
omit  the  derivation  detail  since  it  parallels  that  for  the  one  dimensional  case. 

We  have  the  following  two  remarks: 

(1).  If  we  try  to  do  the  same  estimate  for  a  general  convex  entropy  U(u),  we  can  obtain 
exactly  the  same  cell  entropy  inequality  module  an  interpolation  error  term: 


(2.12) 


-  /  {u(x,t)t  +  /(u(x,<))*)((/'(u(x,f))  -  UU’(u(x,t)))dx 

where  FIv  is  a  projection  into  the  space  V/,  which  interpolates  at  the  two  end  points  of  1:. 
This  motivates  the  following  /i-dependent  modification  to  the  scheme,  which  is  similar  to  the 
“shock  capturing  term”  added  to  the  streamline  diffusion  method  in  Johnson,  Szepessy  and 
Hansbo  [9]  and  in  Cockburn  and  Gremaud  [6]: 


/  ut(x,t)v(x)dz  +  hj+i(t)v(x.  i)-~h-L(t)v(x+l)-[  f(u(x,t))vx{x)dx 
Jl}  J  J  2  J  2  Jlj 

+  Ch  [  \ut(x,t)  + f(u{x,t))x\  ^f^-vx{x)dx  =  0  (2.13) 

Jl,  t)\ 


where  C  is  a  suitable  positive  constant.  This,  together  with  a  bound  which  can  be 
obtained  by  using  nonlinear  limiters  [3],  will  give  us  a  cell  entropy  inequality  for  arbitrary 
convex  entropy  at  least  for  *he  r  =  1  (second  order)  case,  hence  convergence  for  any  noucon- 
vex  flux  f(u)  in  this  case.  For  general  r,  more  ^-dependent  modification  is  needed.  Notice 
that  by  adding  this  modification,  the  formal  order  of  accuracy  of  the  scheme  is  not  changed: 
|ut(x,<)  4-  f(u(x,t))x |  is  just  the  local  truncation  error.  A  similar  argument  as  in  [9]  shows 
that  the  modification  does  not  destroy  convergence  towards  weak  solutions  (conservation). 

(2).  The  so-called  discontinuous  Galerkin  method  can  also  be  recasted  as  a  finite  differ¬ 
ence  scheme  (Hermite  type:  where  one  evolves  both  the  mean  and  the  slope,  maybe  more). 
For  example,  the  second  order  case  is  just  the  following  scheme  for  the  mean  u}(t)  and  the 
slope  Sj(t)  where  u(x,f)  =  Uj(t)  +  Sj(t)(x  —  Xj)  in  cell  1,  (see  [3]): 


2l3>(‘)  =  [/*+}(')  + /<-}(')]  +  2&/  /(MO  +  M'H* -*)))*  (214) 

where 

/,+  *(<)  =  *  (“,(<)  +  ^», (<),“,«(<)  -  (2.15) 
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We  have  thus  in  effect  proved  cell  entropy  inequalities  for  such  high  order  Hermite  type 
finite  difference  schemes  without  using  the  help  of  any  nonlinear  limiting.  Of  course,  to  get 
convergence  one  must  use  the  nonlinear  limiting  to  obtain  L and/or  total  variation  bounds. 


3  Time  Discretization 

We  discretize  (2.1)  in  time  by  the  following  class  of  methods: 

/;,  U  (JA<  U  (X)v(x^dx  +  ^+|u(xJ+i)  ~  h¥\v(*U) 

-  /  f(u"+9(xj)vx(x)dx  =  0  (3.1) 

•'b  / 

where 


un+9(z)  =  (1 -0)un(x)  +  0un+1(*) 


For  6  =  0  this  is  the  Euler  forward  discretization;  for  6  =  1  this  is  Euler  backward;  and 
for  9  =  j  this  is  Crank-Nicolson. 

If  we  take  v(x)  =  un+9(x),  we  obtain  just  as  before 


Thus  a  sufficient  condition  to  get  the  cell  entropy  inequality 


— ("lMPcIx  +  F’*v  -  Fn+?  <  0  (3.6) 

is  just  6  >  i.e.  implicit  schemes  from  Crank-Nicolson  to  Backward  Euler. 

Up  to  now  we  have  not  considered  the  nonlinear  minmod  type  limiters  in  [3]  and  [4]. 
These  limiters  will  render  the  scheme  total  variation  stable  (in  one  space  dimension)  or  L  ^ 
stable  (in  multi  space  dimensions).  We  have  to  make  sure  that  these  limiters  do  not  destroy 
the  cell  entropy  inequality  (3.6).  This  turns  out  to  be  a  simple  issue:  suppose  u"+1(x)  is  the 
solution  obtained  from  the  scheme  (3.1)  without  limiter,  then  projected  to  get  the  solution 
at  time  level  n  +  1  by  some  minmod  type  limiter  un+1(x)  =  Pun+1(x)  (this  is  the  procedure 
adopted  in  [3]  and  [4]).  The  above  derivation  for  the  unlimited  scheme  will  give  (3.6)  for 
un+1 : 


Ji,  2A t  >+?  i-t  -  v 

hence  a  sufficient  condition  to  get  the  cell  entropy  inequality  (3.6)  for  the  limited  (projected) 
solution  un+1  is  to  require  the  projection  P  to  satisfy 


(3.8) 


for  all  polynomial  iu(x)  of  degree  up  to  r  (recall  that  r  +  1  is  the  order  of  the  scheme),  i.e. 
P  does  not  increase  the  L2  norm  of  any  r-th  order  polynomial  in  cell  /_,.  Notice  that  this 
is  the  idea  used  in  [2].  An  easy  way  to  ensure  condition  (3.8)  is  to  write  w(x)  and  Pw(x) 
as  expansions  of  (scaled)  Legendre  polynomials  p;(x)  as  in  [3]:  tu(x)  =  £J_0  aiPi(x)  and 
Pw(x)  =  A  simple  sufficient  condition  to  ensure  (3.8)  is  now  just  |6/|  <  |n/| 

for  all  /,  which  is  easily  checked  to  be  correct  for  the  P1  case  (r  =  1)  in  [3]  and  also  for  the 
general  PT  case  with  similar  minmod  limiters  on  a*  (see  also  [1]  for  such  limiters). 

It  is  also  possible  to  get  cell  entropy  inequalities  for  certain  explicit  time  discretizations. 
For  example,  with  leap-frog  time  discretization: 


8 


I, 


uw+,(*)-un“1(x) 
2A  t 


v(x)dx  +  hj+i_{t)nv(xj+k)  -  hj-}(t)nv(xf_i)  ~  Jj  f{un{x))vx{x)dx  =  0 

(3.9) 


we  can  obtain,  by  taking  v(x)  =  u"(x): 


L 


un+1(x)un(x)  -  un(*)«n~,(x) 

2  At 


dx  +  F"  , 

JT  2 


-FT  i  <  0. 


4  Concluding  Remarks 


(3.10) 


We  have  shown  that  the  discontinuous  Galerkin  method,  or  equivalent  Hermite  type  fini‘e 
difference  method,  which  evolves  the  whole  polynomial  inside  a  cell,  satisfies  a  cell  entropy 
inequality  for  the  square  entropy  in  mu'fi  space  dimensions.  Limiters  to  enforce  total  varia¬ 
tion  or  maximum  norm  stability  can  be  designed  not  to  destroy  this  cell  entropy  inequality. 
Time  discretization  is  also  discussed. 
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